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Matrix model version of AGT conjecture and 
generalized Selberg integrals 

A.Mironov* Al.Morozov^^ and And.Morozov^ 

Abstract 

Operator product expansion (OPE) of two operators in two-dimensional conformal field theory includes a 
sum over Virasoro descendants of other operator with universal coefficients, dictated exclusively by properties 
of the Virasoro algebra and independent of choice of the particular conformal model. In the free field model, 
these coefficients arise only with a special "conservation" relation imposed on the three dimensions of the 
operators involved in OPE. We demonstrate that the coefficients for the three unconstrained dimensions arise 
in the free field formalism when additional Dotsenko-Fateev integrals are inserted between the positions of 
the two original operators in the product. If such coefficients are combined to form an n-point conformal 
block on Riemann sphere, one reproduces the earlier conjectured /3-ensemble representation of conformal 
blocks, thus proving this (matrix model) version of the celebrated AGT relation. The statement can also be 
regarded as a relation between the 3j-symbols of the Virasoro algebra and the slightly generalized Selberg 
integrals ly, associated with arbitrary Young diagrams. The conformal blocks are multilinear combinations 
of such integrals and the remaining part of the original AGT conjecture relates them to the Nekrasov 
functions which have exactly the same structure. 



1 Introduction 

The AGT conjecture [1] unifies and identifies a number of different domains in modern theory, what makes it 
a very interesting and promising subject, attracting a lot of attention [2]- [45]. In its original form the AGT 
conjecture relates the conformal blocks in 2d conformal field theory (CFT) [46, 47] and the Nekrasov functions 
[48], obtained by expansion of the LNS multiple contour integrals [49]. In this form it is now proved only in 
three cases: in the limit of large central charge c [19], when the conformal blocks and the Nekrasov functions 
reduce to (generic) hypergeometric series; in the case of special value of one of the external dimensions [14, 6], 
when they are also hypergeometric series (however, different from the first case); and in the case of a 1-point 
toric function [33], when one can use the powerful Zamolodchikov recurrent relation [50, 22]. Following earlier 
considerations in [14, 11, 29, 30, 31, 32] in [39] a simpler version of the AGT relation was suggested, identifying 
conformal blocks with the Dotsenko-Fateev /3-ensemble integral [51, 52, 58, 11], which can be considered as a 
new avatar of the old proposal in the free field approach to CFT [51, 53] and, at the same time, as a concrete 
application of the more recent theory of Dijkgraaf-Vafa (DV) phases of matrix models [54] . Ref. [39] contains 
absolutely explicit formulas for generic conformal blocks, made from the free field correlators with screening 
integral insertions, analytically continued in the number of screenings. The only problem is that these formulas 
are very tedious to derive and their meaning from the point of view of representation theory of the Virasoro 
algebra, an underlying algebra for the standard construction of the conformal blocks in [46, 47], remains no less 
obscure than in the original AGT relation of [1] . 

In this paper we provide a similar, but conceptually different derivation of the same formulas of [32, 39], 
which involves nothing but the Virasoro representation theory and by now the elementary Selberg integrals 
[55]. Calculations remain tedious but now they are conceptually clear and straightforward. We give only basic 
examples, but a full constructive proof can definitely be worked out in this way, and in this sense one may say 
that the simplified version [39] of the AGT conjecture [1] is now practically established. In fact, the conceptual 
proof is readily available and is given in s. 11.1 in the conclusion. It can be further promoted to a straightforward 
proof of the original AGT conjecture: after the recent progress in [45] there remain just a few combinatorial 
details to fix. In this sense the program to prove the AGT conjecture through the technique of the 

* Lebedev Physics Institute and ITEP, Moscow, Russia; mironov@itep.ru; mironov@lpi.ru 
'fITEP, Moscow, Russia; morozov@itep.ru 

^ITEP, Moscow, Russia and Physics Department, Moscow State University, Moscow, Russia; andrey.morozov@itep.ru 



Dotsenko-Fateev (matrix model like) integrals, which was formulated in [14] and [11], is nearly 
completed. Still it would be interesting to work out some other proof, establishing an explicit relation between 
the Dotsenko-Fateev and LNS integrals, perhaps, making use of the duality between Gaussian and Kontsevich 
models: this is, however, only mentioned in s.11.2 in the Conclusion and remains beyond the scope of the present 
paper. 

In this paper wc consider the triple functions, the coefficients in the operator product expansion involving 
a sum over Virasoro descendants, in the free field formalism. We begin with the case of the single free field, 
i.e. with the pure Virasoro chiral algebra (in the AGT terminology this corresponds to the {7(2)®^""^^ case), 
the extension to k fields, the Wk+i chiral algebra and the U{k + 1) quiver is straightforward. In tcirms of the 
Virasoro primary, elements of the Verma modules look like : L-yVa '■ and the operator product expansion in 
the corresponding chiral algebra is 

: L.rMO) ■■ : L.y^V^M : = ^ A. E ''"''^trn; A.,r. : L.yV^{0) : (1) 

A Y 

Here Va denotes the Virasoro primary with dimension A: i„VA = for n > 0, LqVa = AVa, and L-y = 
...L-njL-nj, denotes the "negative" (raising) Virasoro operator, labeled by the Young diagram Y = {ni > 
712 >. . •}• The coefficients S^^^^ depend on choice of the conformal model. In particular, their values, directly 
provided by the free field formalism below, are usually referred to as the Liouville model structure constants 
[56], we do not consider other choices in this paper. In contrast to 5*^^^^, the coefficients C^'^Yi- A2 Y2 ^'"'-^ 
universal, depend only on the properties of the Virasoro algebra, and these are the quantities we are going to 
investigate. Moreover, we further restrict our consideration to the case of 12 = 0, this is enough to reproduce 
the spherical 4-point conformal blocks, studied in [1, 35, 39]. 

The coefficients C can be straightforwardly found by standard CFT methods, see [5] for a detailed review. 
Coming back to free fields, they provide an alternative derivation, somewhat simpler and more transparent, see 
[5] and s.2. The only problem is that in the free field model there is a "conservation law": the primaries are 
represented as Va = e"* with^ 

A = a{Q-a), c=l-6Q^, (2) 

and in the sum at the r.h.s. of (1) 

a = ai + a2 (3) 

Thus, only a restricted set of the triple functions C, namely, C'^^ y-^.^'^^ y^ can be defined in this model (from 
now on we label these functions with a- rather than A-parameters). It is a long-standing problem in CFT, how 
the free field formalism can be used to obtain arbitrary C^'^y^- a2 Y2 '^ith a ai + a2. The results of [32, 39] 
imply that the operator product 



(4) 



^Iii terms of paramctrizatioii from the Appendix of [6], normalization conventions for free fields in the present paper are as 
follows; fei = 1/2, fc2 = 2, p = 1/A = 1, where 

</-(2)<^(0) = k-z log 2, 
2,K2 

c = 1 - 12k'ik2Q\ 
A(ya) = ^a(a-Q), = 2fci Q/p 

Then 

■■ L_2Vc : =: W)^ + (fciQ +pa)a2</.j Vc, : 
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where b is the Dotsenko-Fateev screening charge, i.e. Q = b — 1/b, has exactly the same expansion coefficients 
as OPE, 



(5) 



and we demonstrate below that this is indeed the case. Eq.(5) is the main claim of the present paper, supported 
by a number of examples. In other words, the r.h.s. of eq.(4) is identically the same as (1) provided the structure 
constants C and C are related by a change of variables (2) and additionally 

a = ai +a2 + bN (6) 

Thus, (4) resolves the above mentioned problem in the sense of analytical continuation: the coefficients 
C^'^Yi- a2 Y2 rational functions of a and they are fully defined by their values at discrete points a = 
ai + a2 + bN. 

If the original two fields are primaries, Yi =¥2 = 9, then eq.(4) is derived in three steps. 

A) First, one uses the basic free field relation, 

N ( N N 1 

2=1 I i<j i=\ J 

and then expands the exponential in powers of q and Zi : 

. ^oi^m+cc^<P{i)+bT.i<P{^i) . = ^ q\^\-\^'\HY'YZ^' : L_Ye''"^+"'^+^^^^^°^ : (8) 

Y,Y' 

Here = J^,- z"' for a Young diagram Y' = {ni > n2 > . . .} and the sum goes over all pairs of Young diagrams 
with \Y\ > \Y'\. At this step, one evaluates the z-independent coefficients Hyy', as functions of ai,a2 and N. 

B) Next, one takes the integrals over Zi, 

iy' = n /' dz. I n(-^ - ^^r^ n - -o^^-^ | (9) 

i=l'^° [ i<j i=l J 

For the single-line Young diagrams Y' = [!"] these are the well known Selberg integrals, which generalize the 
Euler S-function and are equal (after the standard analytical continuation from integer powers in the integrand) 
to the ratio of Gamma- factors. For generic Y', the integrals generalize the Selberg integrals producing extra non- 
factorizable polynomial factors, which can be explicitly evaluated. Being polynomial, they do not complicate 
the analytical continuation. 

C) Combining the results of steps A and B, one gets the structure 



^aua2 = X] HyY'Iy' 
\Y'\<\Y\ 



(10) 

9=1 



in the form of a finite sums over Young diagrams. 

D) The last step is to compare the C with the known expressions for the conformal theory structure constants 
C (the 3j-symbols of the Virasoro algebra), transformed with the help of (2). 

For non-trivial Yi and Y2 the calculation goes the same way, with additional powers of z^^ and {zi — q)~^ 
emerging in the integrand. 



In this letter we provide in full detail a sample calculation of this kind for the two simplest cases of 
{Yi,Y2,Y} = {[0], [0], [1]}, {[0], [0], [2]} and {[0], [0], [11]}. It is enough to demonstrate the principle and can be 
straightforwardly computerized to provide more examples. There is small doubt that all such examples would 
confirm the relation, which at the moment looks like a non-trivial statement, identifying the 3j-symbols of the 
Virasoro algebra with linear combinations of the generalized Selberg integrals Iy- 
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2 The free field formulas for C^'+'^' 

Evaluation of the operator product coefficients in the free field model is considered in detail in [5] . The simplest 
example is: 



_ ^a\a.2 



Now, using 



from 



1 + 



,(ai+a2)0(O) 



(11) 



(12) 



one obtains 



c: 



'cei+a2, L-i 



a2 



^aia2 



^OC\OC2 



Aal + Q!2 (2a2(5 + 4q:iQ!2 - 1) 
2 (^4 (ai + a2f + 2Q (ai + a2) - l) (ai + 0:2) 

2aia2 



4 (ai + Q!2) + 2g (ai + 0:2) - 1 



(13) 



3 The CFT formulas for C"^ 

According to [47, eq.(5.16)], for three generic dimensions one has, instead of (13), 



^0.10.2 



A2 + A - Ai 
2A 



Similarly, 



a,Lii _ 8A3 + cA^ + 16A2A2 + 2AA2C - leA^Ai - 2AAic - 4A2 + cA + SAjA 
~ 4A(16A2-10A + c + 2cA) ^ 

A^c - I6AA1A2 - 2Ai A2C - I6AA2 + A2C + 8Af A + Afc + 4AAi - Aic 



+ 



4A(16A2 - lOA + c + 2cA) 



(14) 



(15) 



L_2 _ A2 + 2AA2 + 2AAi - A - 3A2 + 6A1A2 + A2 - 3A2 + Ai 



16A2 - lOA + c + 2cA 

and so on. 

Only in the case of the conservation law condition, 

a = ai + a2 

these expressions are reproduced by the free-field formula (13). 



(16) 



(17) 
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4 Operator product with screening insertions 

In order to relax the 17(1) conservation law (17) we insert the Dotsenko-Fateev screening charges into the l.h.s. 
of (11). Then, in the integrand one has 

N N N 

: e"'"^^°^ : : e"^'^*'^ : : e''*''^''*^ : = g2"i"2 ^^(0^ - Zj)'^^'' JJ-z^''"'(g - Zi)^''"^ : e"i'^W+"2'^(9)+''^< "^(^'^ : 

1 i<cj i=l 

N N N 

i<ij i—1 i—1 

( .V^ V (^^(0))' ( 2 2^ d^Q) 



3\ 930(0) 



3! 



+ ... e 



(ai+a2+6Af)0(O) . _ 



_ J2oiioi2 Y\ I ^ ^ ^26^ TT ^25ai / ^ \2ba2 ■ / i i '^2g + 2^^=1 ^' r , /i q\ 



(4(ai + a2 + 6Af) + 2g) (aag + 6 Eti - («2'7' + & Eti 
^ 2(ai+a2 + 6jV)(4(ai+a2 + bA^)2 + 2Q(ai + a2 + ^'iV)-l) ^"'^ ^^^^ 

(ai + a2 + 67V) (a2g' + ^Eti - («2« + &Er=i ^if 

4(ai + a2 + 6A^)2 + 2Q(ai + as + &A^) - 1 ^"^"^ 



AL-li + 2(q;i + a2 + bN)BL^iL^2 + 2(ai + 02 + hN)CL^^ 



(20) 



6(ai + a2 + 6Ar)(4(ai + 02 + 6iV)^ + 2Q{ai + 02 + &iV) - 1)((q;i + 02 + 6iV)^ + Q(ai + Q2 + hN) - 1) 

_|_ -J g(ai+a2+feA')0(O) . 

where 

A = i^2q^ + ^E ^i* j ~ ("1 + 0:2 + + Q) ^0:25 + ^Yl^^ ^"25^ + j + 

+ (4(ai + a2 + 6A^)^ + 6Q(ai + a2 + 6A^) + 2Q2(q,^ + a2 + feiV) - 2) |^a2g + ''E-^ij 

B = -4(ai + a2 + 67V)2 ^a2g^ + & ^ ^'"^ j + 
+4(ai + a2 + 6Af)2((ai +a2 +6iV) + g) ^025 + 6^-^*^ ^"25^ + 6 E -^i^ j + (21) 



+4 



(1 - 3(ai + a2 + 6iV)2 - 3Q(q!i + a2 + hN)^ (^2q + ^^1^^ 



C = 2(4(ai + a2 + bNf + 2Q{ai + a2+hN)-l) {{ai + + bNf (^2q^ + h ^ z| j + 
Underlined by one, two and three lines are the contributions at levels one, two and three respectively. 

5 Structure constants C^^^^"^^^^ from Selberg integrals. Level one 



Now, as suggested in [32, 39], we take integrals over Zi along an open contour which connects positions of the 

two original operators. In order to perform integration at level one, i.e. to integrate the first line in OPE, (18), 
one needs the integrals which are given by the now standard formulas from ref.[55] (see also the Appendix in 
the present paper): 



.-1 "'0 



jv-1 rfl + 26ai + jbAvU + 2ba2+]bAT(l + (j + l)bA 
j}o r (2 + 26ai + 26q2 + {j + N - 1)62^ r (l + 

and 

N „a N N 



j=l "'O j<7 i=l 



TV „1 JV N 

^^N+1){iWN) . ^2(ai+JVb)(a2+W6) fT / - r.^26«2 TTCr. _ ^.^26= /'^^ j_ ^ \ (70)&(71) 



^(JV+l)(l+6''iV) . ^2{ai+Nb){a2+Nb) 

3=0 



n / dz,zf'"^^{l-z,r"^ll{z,-zjr'(a2 + bY,z?) 

i=l i<j i=l 

r(l + 26ai +i62)r(l + 26a2 +i62)r(l + (j + 1)52) 

r(^2 + 26ai + 25a2 + (iV - 1 + jW^jV^l + b^'j 



1 + 2bai + (N - l)b^ \ 
"-+^^ 2 + 26a. + 26a2 + 2(Ar-l)62 j ^''^ 

One can now extract the structure constants from the integrals of (18). We do it first for A'' = 1 where formulas 
are just a little simpler, and then for arbitrary N. 

• For N = 1 these integrals are just the Euler i?-functions: 

<i>=rdz z^^-^ (1 - zr-^ = ^^ix '':f^' ^ 

Jo r(2 + 26ai + 26a2) 

^ ^, ^ 25awi ^26a2/ r(l + 26ai)r(l + 2ba2) / ^ , l + 26a^\ 

< a2 + bz > = dz z ^ (1 - ^) ^ (a2 + ^-2) = ' — — ^ 02 + 0- — — — — (24) 

Jo ^ ^ — T{2 + 2bai+2ba2) \ 2 + 2bai + 2ba2j ^ ' 

According to (4) and (18), the first of these formulas defines 



cai+a2+6 (2J r(l + 26ai)r(l + 2ba2) 
r(2 + 2bax + 26a2) ' 



(25) 



while the second one is proportional to the product S'^\a'^'^^^Ca\at'^^^'^ ^ ■ Therefore, Cala"^^'''^ ^ is given 
by the ratio of the two integrals (up to an additional factor): 

g,«i+a2+b.L_i _ <a2 + bz> 1 (24) / 1 + 26ai \ 1 

<1> ai+a2 + b V 2 + 2bai + 2ba2 J ai+a2 + b ^ ' 
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At the same time, from (14) in this case one has, taking into account that Q = b— 1/6: 

^a,+a2+b,L_i ^ A2 + A - Ai (J a2{a2 -Q) + {ai +a2 + b)iai +a2 + b-Q)- aijai - Q) ^ 

2A ~ 2{ai + a2 + b){ai+a2 + b-Q) 

2a2{ai + a2) + 2bai+2a2/b+l ( , 1 + 26ai \ 1 

— ' a2 + 



2(ai + a2 + &)(«! + a2 + 1/6) V ^ + 26ai + 26a2 / ai + a2 + 6 

(24) <a2^bz> 1 (26) ^„,+a,+6,I,_i 

< 1 > ai + a2 + 6 



= C'-^r+''''-^ (27) 



Similarly, for arbitrary N: 



^a,+a2+bN,L_r ^ Ag + A - Aj ^ a2(a2 -Q) + jai + Q!2 + bN){ai +a2 + bN -Q)- aijai - Q) ^ 

2A 2{ai+a2 + bN){ai+a2 + bN -Q) 

+ h^r ^ + '^ba^ + {N-l)b^ \ 1 ^ <a2 + bEiZi> 1 ^ ^„,+„,+6iv,L_. .^sl 

^ 2 + 2bai + 2ba2 + 2{N - l)b^ J ai+a2 + bN < 1 > ai+a2 + bN "'"^ ^ ' 



and 

^irfl + 26ai+j62yfl + 26Q2+j6') Hi + 762^ 
/=o r (2 + 26ai + 26a2 + iN-l+ i)^^) fj, r(l + 62) 

Formulas for the structure constants C are rational, therefore, they can be straightforwardly analytically con- 
tinued in N to arbitrary values of a = ai + 02 + 6A''. The analytical continuation of the above expression for 

5 is somewhat more ambiguous (and, anyway, there is nothing to compare them with, since the coefficients S 
generally do not factorize into holomorphic and anti-holomorphic parts). 

6 Level two 

At level two, one needs to integrate (20). The ordinary Selberg integrals (70) and (71) are not sufficient for this 
purpose, one needs also the generalized one (73) from the Appendix. Then, the two integrals that one needs in 
(20) turn out to be 

N .q N N 

«^--n / dz.zr''^iq-z,r'^^l[{z,-z,r\a2q' + b^z^) = 

i=l i<j i=l 

AT 1 N N 



= ^(iV+l)(l+6^Af) + l . q2{ai+Nb){a2+Nb) "Q 



n/ dz,zr'{l-z,r'^^ll{z.-z,r\a2 + bY^z!) 

i=l "'O i<j i=l 

a2 + A^6 X (30) 



<3 

N-} r(l + 26ai + j62)r(l + 26a2 +i62)r(l + (j + 1)62 



,=0 r(2 + 26ai + 26a2 + {N - 1 + j)62)r(l + 62 
(4a2b2 + 40,^0,26^ + 6aib^N - Saib^ + «a^a-2.N - 4&3a2 + 8aib + 4a2b + 4 + Sfc^ATZ _ y{,4^ ^ 4{,4 _,_ 752 _ 962)(2ai6 + \?N -\? -V\) 



2(2ai6 + 2026 + 2b2W -362 + 2)(2ai6 + 2a26 + 2b2Ar - 26^ + 3)(ai6 + 026 + b^JV - b^ + 1) 

and 



N .q N N 

n / dz,zr''-^{q-zr''"'l[{z,-zjr"'\a2q + bY,z,) 

i=l "'O i<j i=l 

N N , ^ . 2 



^(JV+1)(1+6=JV) + 1 . g2(a,+iV6)(a.+iV6) "Q f dz^zf'^- {I - Zif>"^^ f[{Zi - Zjf'^ (a2 + & "^''^ 

1=1''^ i<j i=l 



N-1 



t(i + 2bai + jbArll + 2ba2 + jbAvfl + {j + 1)62) 

^(JV+l)(l+6=JV)+l . ^2(ai+JV5)(a2+]Vfe) TT _V / V /V ^ ' ) (3^) 

/=o r (2 + 26ai + 26^2 + (A^ - 1 + i)62) F (l + 62) 
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2 nAT u 1 + 2feai + (Af - 1)62 (2ai&+62iv-(,2_,_i)(2ai6 + 62^r_262 + i) 



ai + 2Na2b ' ^— — ^ — + N{N - 1)6^ 



2 + 2bai + 2ba2 + 2{N - l)b^ ^ ' 2{a\b + aib + b"^ N - b"^ + \){2a]_b ^ 2a2b + 2b'^ N - Zb"^ + 2) 

(4Q;f&2 ^40,^0,2^2 + 6a ib^ AT - Saife^ +^b^a2N - 4b^a2 + 8aib + 4a2b + 4 + 3b'^ - 7b^ N + 4b^ + Ib'^ N - 9b'^){2aib + b'^N - b'^ + 1)' 
■ 2(2ai6 + 2a26+ 262JV -362 + 2)(2aib + 2a2b + 2b'^N - 262 + 3)(ai6 + 026 + 62Ar - 62 + 1) ^ 



7 Higher levels 

The detailed explicit formulas are quite lengthy already at level two. Writing them down for higher levels is 
simply impossible: they take several pages. However, in every particular case eq.(4) can be easily validated by 
simple computer calculations, provided one knows the following set of matrices. 

• The action of Virasoro generators on the free field primaries, 

L-ye""^ = ^YY'{a) : J^' e°"^ : (32) 

\Y'\ = \Y\ 

where = d"^ (p d"'^ (/) . . ., produces the matrix Cyv (see s.2). 

• The expansion 

Y 

gives the vector Sy actually expressed through the Schur polynomials. Up to level 3 this vector, and also 

£:^= ^ L^y'Sy, (34) 

\Y'\ = \Y\ 

with being the inverse of Cyy', are explicitly given in eqs.(18)-(21). 

• If Yi and Y2 in (4) are non-trivial, then one actually needs a more sophisticated triple- vertex 
describing the expansion 



: J^^e«^^(°) : : J^=e«=^(«) : : e''^^^^^^) := ^ £^^^=(g, 5) : jy(o)e(«^+«=+''^)^W : (35) 

Y 

In this case, one has to consider also the quantity 

SnY.;y= ^ Cy^^,Cy^y,S^)''iC^^' (36) 

YI,Y^,Y' 

• These £y are actually functions of {zi}, i.e. have the form 

£y = Y^ Syy z^' (37) 

Y' 

with one extra index Y'. This time the sizes of Young diagrams can be different, only |y| < \Y\. 

• Integration over z converts £ into 

<£y > =Y^£yy,I^ (38) 

Y' 

where = < z^ > are generalized Selberg integrals, described in the Appendix below. In fact, as emphasized 
in [45], they are expressed through the simpler quantities, the averages of Jack polynomials with the help of 
one more matrix, 

J2 V'''"<Py,> (39) 
\Y'\=\Y\ 

inverse to the matrix of expansion of the Jack polynomials into monomials, 

Py =Y^VyY,Z^' (40) 
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Putting all the things together, one obtains for the Dotsenko-Fateev representation of the conformal triple 
function: 



Gy,y, = CY,Yi{ai)CY,Yi{a2)C^''' {ai + as + hN) < £yV' >, 



■^Y'Y' 



>, 



> 



(41) 



Summation over repeated indices is implied. 

This should be compared with the usual CFT expression 



YY' 



(42) 



Details of this calculation are described in [5]. 
An explicit check of the relation 



Y1Y2 



c 



Y1Y2 



(43) 



at levels 1 and 2 and for Yi = I2 = is described above in ss.5-6. Since all the matrices are explicitly presented 
there also for the case of level 3, it is a trivial computer exercise to make a check also at this level, and, of 
course, it also confirms relation (43). 

To check it at other levels for Yi = Y2 = $ one needs to know four matrices, C, £, V, Q, and two vectors, 
< > and Ty. When Yi and I2 are non-trivial, £ and T acquire an additional pair of indices, Yi, Fs- These 
entries belong to different sciences: C and £ to the free field calculus, Q and f to CFT, V to the theory of 
orthogonal polynomials, < Py > to the theory of Selberg integrals^. Eq.(43), the weak (matrix model) form 
of the AGT conjecture establishes a concrete relation between the seemingly unrelated quantities from these 
different subjects. 



8 Virasoro intertwiners 

Expansion rule (1) implies that the structure constants C are the components of the Virasoro intertwining 
operator between Verma modules Ai, A2 and A. Comultiplication in the Virasoro algebra is somewhat non- 
trivial [57]: 



A(L_i) = L_i 7 -I- 7 (g) L_i, 
A(L_o) = Lo^I + I^Lo + qI^ L-i, 
A(Li) = Li (g) 7 -h 7 (g) L_i -I- 2g 7 (g) Lo + 7 L_i, 

(44) 

As the simplest example, this means that Lq acts on the operator product expansion of two primaries as 

Lo(VAMVA2{q)) = (^Ai+A2 + g^^ (Fa, (0)FA,(g)) (45) 
This is, of course, in a perfect agreement with (1): 

VAMVA2{q) = 5^gl^lC^'^AiA2i_yyA(0) (46) 

Y 

On one hand, (Ai + A2 -h 5^) q^-^^-'^^+l^l = (A + \Y\)q'^-'^^-'^^+\^\ , on the other hand, LoL_yVa = 
(A + |y|)L_y Va- It is instructive to see how this works also for representation (4). Then, at the l.h.s., one has 
a product of many operators, Vai {0)VA2iq) {Jq Vi{z)dz^ , and the multiple comultiplication (it is associative) 
now acts as 

Lo^I(E)I®^ + I(E){Lo + qL_i) 7®^ + 7 ® 7 ® ((Lq + ^i^-i) ® 7®(^-i) + . . . + ® (Lq + z^L.^j) (47) 

^Note that the Selberg integrals (see the Appendix) produce in the denominators the products automatically presenting the 
decomposition of the Kac determinants in terms of a- variables, [51, 47]. 
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Since Vi — : e'"^ : has unit dimension, {Lo + zL^i)Vi{z) = ^(zVi(z)). Our definition of Selberg integrals is the 
analytical continuation from the points where all 2aj6 and are positive integers, therefore, all ^-derivatives 
are always integrated to zero as if one uses the closed contours, so that one actually gets 



Vi{z)dz] 



(48) 



Note that the (j'-derivative does not act on the upper limit of the integrals. The multiple integral at the r.h.s. 
of this formula depends on q through the factors g2aia2 ^^^^ Yl^^i{q — -^i)^"^'' in the integrand. The action 
of the logarithmic g-derivative on the first factor gives the factor 2aia2, while the action on the others gives 
— 2Q:2&X]j=i When N = 1 this simply means that the Selberg (Euler) integral has its argument c = 2a2& 

shifted by —1: 

2a2bq [" z^'^'^l - zf'^^^-'^dz = (2ai + 2a2& + 1) C z'^"^''{l - zf^^^dz (49) 
Jo Jo 

For > 1 this is a similar, but a little more complicated exercise (see eq.(61) below for a similar evaluation of 

X^i il))' which gives 

2q!2&/ y \ = N(2aib + 2a2h+(N -1)1"^ + 1] (50) 

Substituting all this together with (4) and A, = ai(a, — 6 + 1/6) into the r.h.s. of (48), one obtains for the 
coefficient in front of Va 

al + al- (ai + a2){b - 1/6) + 2aia2 + N{2aib + 2a2b +{N - 1)6^ + l) = 

= (ai +a2 + bN){ai +a2 + bN-b+ 1/6) = A, (51) 

as needed. In a similar way, one can act with any other Lk on (4) and check in detail that it is indeed consistent 
with the comultiplication rule (44). This, of course, follows from the general argument of [51], since the screening 
insertion is an integral of the dimension one operator, and, once again, our definition of Selberg integrals actually 
allows one to consider the integration contour as closed. As we demonstrated in this section, an explicit check 
confirms this general claim. 



9 From OPE to conformal blocks 

Eq.(4) is very well suited for constructing arbitrary conformal blocks. If we denote the operator product 
expansion at the l.h.s. of Eq.(4) through Vi(0) *n V2{q), then conformal block is the value of a linear form on 
an ordered product, for example, 

(((Vi(a;i) *N,2 V2ix2)j *iV(i2)3 Vsixsfj *n^^^^^3^^ V^ixi) * . . (52) 

for Fig.l or 

*Ni2 ^2(2:2)) *JV(i2)(34) (^3(2:3) *Ar34 V4{Xi)^^ *JV((12)(34))5 V5{X5)*...) (53) 

for Fig. 2, analytically continued in all the A'^-variables, which are in this way converted into arbitrary interme- 
diate dimensions. 

The linear form here is defined by the usual rule 

(L_y^e«^(-)) ~(5y,05„,Q (54) 

Note that the product *jv is defined in (4) asymmetrically: the result is an operator at point 0, i.e. at the 
position of the first entry of the product. This makes *jv non-associative: 

(Vi{xi) *N V2{X2)) *M V^ixs) = Vi{xi)V2{x2)Vs{x3) (^J^ ' : e"* :^ (^^ ' : e'^ ^ (55) 
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xi,a\ ai2 = ai + a2 + bNi2 a(i2)3 = "12 + as + ^^(12)3 "((12)3)4 = "(12)3 + "4 + 5iV((i2)3)4 

Figure 1: A comb-like conformal block from [1, 35], for which the AGT relation is known in the case oi = xi ^ xq ^ X3 ^ X4 ^ 

. . . Shown are the o-parameters, the dimensions arc equal to A = a{a — b + 1/6). The intermediate dimensions are parameterized 
by the Af-variables, after analytical continuation in TV they take arbitrary (continuum) values. 



X3,a3 



Xi, ai 





X2,a2 ^ 


/ X5,a5 








"(34) = as + a4 + 67V(s4) 




xi,ai 


Q:i2 = a\+ a2 + hN\2 


a(12)(34) = ai2 + a34 + &-/V(12)(34) 



Figure 2: A star-like conformal block, for which the AGT relation (the corresponding set of Nekrasov function) is yet unknown. 
The "matrix-model" or Dotsenko-Fateev representation, a weak form of the AGT relation, is readily available and provided by 
eqs.(53) and (4) of the present paper. The brackets in subscripts of the iV- variables define the integration segments of the screening 
insertions. Prom the point of view of the conformal blocks, they refer to a particular corner of the a;-moduli space, see [35]. 



while 



Vr{xr) *N (V2{X2) *M Vsixs)) = Vi{xi)V2{X2)Vs{X3) 



: e 



N 



: e 



M 



(56) 



and the difference is in the integration segments in the last items. Thus, the brackets are essential in the 
above expressions for the conformal blocks. In practice, CFT calculations are determined by the ordering of 

.T-arguments: = ^ a::2 <C 2:3 <C 2:4 <C . . . in (52) and = .ti ^ X2 <C <C .x.,, X4 — x^ ^ x^ in (53) and, 
hence, these two cases correspond to different regions of the values of variables in the conformal block. 



10 Towards a proof of the AGT conjecture 

As explained in [32, 39] and further developed in [45], representations like (52) can be directly used to prove 
the original AGT conjecture. For example, for the 4-point conformal block (52) implies that 



^^12 pq " (12)3 p\ 

/s=l"'° j=l i<j k<l i,k 

where q = ^^^^ , xi = 0, X2 = q. xs = 1, Ni = N12, N2 = iVn2l3 and a = 2bai, B = b"^ , c = 2bo.3, 7 = 2ba2. In 
order to make use of the Selberg integrals from the Appendix, which are all along the segment [0,1], we rescale 
Vk QVk and expand in powers of q: 



-^(12)3 \ 



^12 /.I 



n / zr^^'^'^^i-z.rdz.xiz^^ n/ ^^(i-^/^r^y^n^^f 

i=l •'^ i<j fe=l'^" k<l 



fe=i 



(58) 



This nice exponential formula first appeared in [45] . 
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It can be now expanded in the Schur/Jack polynomials so that the result is a bilinear combination of the 
generalized Selberg integrals over the z and y variables. The integrals are labeled by Young diagrams (see the 
Appendix below), thus, one naturally obtains a bilinear expansion in Young diagrams. In the formulation of 
[45], the AGT conjecture is now reduced to the claim that there arc two different expansions of this type: in 
triple vertices and in the Nekrasov functions. Denoting independent averaging over the z and y variables by 
the double angle brackets, we have [45] : 



T.Y,,Y.nYl)Q-\Y^,Y2)nY2) 

• Despite this is already done in [32, 39, 45], for the sake of completeness we explicitly illustrate the situation 
at the first level of the g-expansion. 

At this level, one needs just two explicit expressions for the Selberg integrals from the Appendix: 

^ iVi2J[l] a+(jV,,-l)/3 + l 



and 



^^(12)3/[-l] %2)3/a'-l[l^(->^-^] a' + C+(A^(i2)3-l)^ + l 

^-zj = /[O] = JM = a' 



Here a' = a + 7 + 2Nx^ and 



1 1 \ (A + A2 - Ai)(A + A3 - A4) , , 

a2a3 + 2/3X:-E2/^+^E^+'=E2/^) =^ ^ ^^^^ 

\ i k i k I 



with [39] 



A = a'(a + 7 + (2iVi - 2)p + 2)/2/3 = (a + 7 + 27Vi/3)(a + 7 + 27Vi/3 + 2 - 2/3)/2/3, 

Ai = (a)(a + 2-2/3)/2/3, 
A2 = (7)(7 + 2-2^)/2/3, 
A3 = (c)(c + 2 -2/3)/2/3, 
A4 = (a + c + 7 + 2/3(7Vi + N2)){a + c + 7 + 2/3(A^i + A^a) + 2 - 2;3)/2/3 (63) 

thus reproducing the expression for the g-linear contribution to the conformal block [46, 47, 4] . 

• In general the q^-term of the B-expansion contains bilinear combinations of the integrals /atJI'] and 
/jVa [—Y'] with \Y'\ < m. Generalization to the multi-point conformal blocks with the multi-linear expansion 
in Young diagrams is also straightforward. 



11 Conclusion 

In this paper we justified the claim that the coefficients of the operator product expansions in arbitrary conformal 
theory are fully controlled by the free field model, provided one allows insertions of the Dotsenko-Fateev screening 
operators between the points, where the original operators are located, and analytical continuation in the number 
of these insertions. The well-known complexity of the operator expansion coefficients appears related to that of 
the generalized Selberg integrals, which are defined for arbitrary Young diagrams, but contain non-trivial non- 
factorizable polynomial factors when the diagrams are different from [1"]. Since the single line diagrams [1"] 
are associated with the hypergeometric series [14], one may say that the non-triviality of the Selberg integrals 
for other diagrams is responsible for the deviation of the conformal blocks from the hypergeometric functions 
and, thus, it is what requires the generic Nekrasov functions to appear in description of the conformal blocks. 
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11.1 The proof of the matrix- model version of AGT conjecture 



Despite the present paper does not contain a full constructive proof, hopefully, it provides a conceptually clear 
explanation of the week form of the AGT conjecture [32, 39], identifying the conformal block with the analytically 
continued matrix model partition function in the DV phase [11, 30, 29, 31]. Moreover, for the 3-point functions 
this identification can be implicitly (not constructively) proved with the following chain of arguments: 

• The structure constants C^'^Yi-a2 Y2 (-'■)' ^-^^ components of the intertwining operator, are unambiguously 
defined by representation theory of the Virasoro algebra. 

• Free field + DF induced structure constants C'^'^Yi-a2 Y2 C'^) components of the Virasoro 
intertwining operator, but for a triple of concrete and explicitly realized Verma modules. Instead they are 
defined only for discrete values of a = ai + 0:2 + bN. 

• The both C and C arc rational functions of their arguments a and a. The rational analytical continuation 
in a of the function C is unique, therefore, such an analytical continuation coincides with C: 



(64) 



C. 



ai,Yi;a2,Y2 



rcx,Y 

'-^ai,Yi;a2,Y2 



Note that with this technique one obtains "matrix-model" representations for arbitrary conformal blocks, 

not only for (52), but also for (53). At the same time, the literal AGT relations are currently applicable only for 
the case of (52), their generalization (an extension of the set of the Nekrasov functions) to arbitrary conformal 
blocks remains to be found. 

11.2 Towards a proof of the remaining pcirt of the AGT conjecture 

After the matrix model version of the AGT conjecture is proved, the original AGT conjecture is reduced to 
an exercise, outlined above in section 10. There are still some combinatorial identities to be proved in this 
direction, but this seems rather straightforward. F complete proof of the AGT conjecture on this track 
is now clearly within reach. 

Much more interesting would be to prove the AGT conjecture differently, without any direct use of the 
Ncikrasov functions. Given the result of the present paper, it turns into a puzzling observation that the two 
kinds of integrals are identical: the matrix model (Dotsenko-Fateev) integrals like (57) and the LNS ones [49] 
like 



k k 



k k\Wj l^^y 27rz(x2-(ai2 -e/2)2)((:.„ + e)2-(ai2- 6/2)2)^^x2, -e?)(x2,-e2) ^ > 

where ai2 = ai + a2 + bNi2, e = ei + £2, ei = b, £2 = —1/6 and the four parameters m^, i — 1...4 are linear 
combinations of ai.2,3 and 0:4 = Qi2 + ck3 + &-^(i2)3 [4]- Thus, the positions of poles in the LNS integral (dictated 
by OL\2 and, hence, by N\2) become a number of integrations in the DF case, while the number of integrations 
in the LNS integral depends on the degree of expansion into g, i.e. on the level of expansion of the conformal 
block into descendant contributions. 

Thus, the AGT acquires form of a duality relation, where the number of integrations on one side is a 
parameter in the integrand at the other side and vice versa. This type of duality may seem mysterious, but 
it is well-known in the theory of matrix models [58]. The simplest example is provided by conversion of the 
Gaussian model into the Kontsevich type model [59] : 



/ dMe-^^^'+^^*'=^^' ~ / dX(detX)^e- 

JnxN Jnxn 



UrX^-itrAX 



(66) 



where Tr , Det and tr , det denote the traces and determinants of the NxN and nxn matrices, and tk = -^tr A"*^ . 
Both models arc of the eigenvalue type and clearly the number N of integrations over the eigenvalues at the l.h.s. 
appears just as a parameter at the r.h.s., where the number of integrations is a fully independent parameter n. 
This identity, (66) can also be rewritten as 



/ rfM Det(/(g)/-M(g)A-^) e""^^' ~ / dXdet(l-i^) 

J NxN ^ ' Jnxn \ 



N 



e-^f^^' (67) 
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which, if expanded in powers of A ^, becomes an identity for the Gaussian correlators: 

11 1 / 1\^ 

1- -tr^ «TrM2»jv+ - (tr-j « (Tr »^ + . . . = 

\ In 

Since << MijMki » = SjkSu and < XabXcd > = S^c^ad, the two sides of the cquahty coincide, but, as usual 
for dualities, the second term at the l.h.s. is equal to the third term at the r.h.s. and vice versa. 

Eq.(66) can be used (at /3 = 1) to further transform the multiple integrals (57) so that N becomes a 
parameter in the integrand. Indeed, such integrals arise from the l.h.s. of (66) for tk = X^a which for 

the integer values of correspond to A matrices of the block form A = ^^^qalixa at the r.h.s. of (66). The 
number of integrations at the r.h.s. is then equal to n = X^aMa- It is natural to conclude that the AGT 
relation between (57) and (65) is a further generalization of the duality relation (66) continued to 
the /3-ensembles (to (3 1) and to non-integer values of n. Details of this analysis will be presented elsewhere. 



11.3 Extension from Virasoro to W 

Another mystery associated with our result in this paper concerns extension to the case of several free fields. 

In conformal theory, this corresponds to switching from the Virasoro to W chiral algebras. The problem is 
that representation theory of the W algebras is not sufficient to specify unambiguously arbitrary conformal 
blocks. Additional constraints should therefore be imposed by brute force. It is, however, unclear what are 
the parallel restrictions in the free field formalism and its Dotsenko-Fateev extension described in the present 
paper, which seems easily generalizable to an arbitrary number of free fields. This subject also remains open 
for future investigation. 
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Appendix. Selberg integrals and their generalizations 



The Selberg integrals 

N N 



lY = lir dz, J f[{zi - zjf^ n zt{q - z,r 

j=l"'0 I i<j i=l 



(69) 



with = z"^^ Z2^ ... for y = {rii > n2 > . . .} arc direct generalizations of the Eulcr Beta- function, also 
represented as products of elementary Gamma-function factors. The Selberg integrals are naturally labeled 
by Young diagrams Y, and well-known are only integrals for the single line diagrams [1"]. For more compli- 
cated diagrams, the integrals contain additional polynomial factors, which arc not further factorized into linear 
expressions. However, they are needed for comparison with DF 3-point functions in this paper. 

• From [55] one knows that 

fJl r(/3 + 1) fj-^ r(a + c+{N-l+j)P + 2) 

III] = - + + ^ m (71) 

^ ^ a + c+{2N -2)p + 2 ^ ^ ^ ' 

and, in general, 

<-) 

• If the Young diagram Y contains fc > 1 lines, then [55] is not sufficient, and actually /[y] acquires additional 
factors, which are polynomials of degree 2k — 2. In particular, 

^ + ac -h (3^- - 4)a/3 + 2[N - l)/3c + 4a + 2c + 4 + (TV - l)(3jV - 4)/3^ + (7A^ - 9)/3 ^^^^ ^ 
(a + c + (27V - 3)/3 + 2) (a + c + (27V - 2)/3 + 3) 

fa + (TV - l)/3 + 1) ia^ -h ac + (37V - 4)a/3 -h 2(7V - l)/3c -h (7V - 1)(37V - 4)^^ + 4^ + 2c (77V - 9)/3 + a\ 

= A Ls L /[o] (73) 

(a -I- c (27V - 3)/3 2^ (^a c -|- (27V - 2)/3 ^ 2^{a ^ c ^ {2N - 2)/3 

(a + (7V - 2)/3 + 1) (a + (7V - l)/3 + l) /[O] 
(a c -I- (27V - 4)/3 + 2^ (a + c (27V - 3)/3 2^ (a c -|- (27V - 2)/3 + 2^ (a + c (27V - 2)/3 -h 3^ 

• {a^ + ac + (37V - 5)a^ + (27V - 3)/3c + 3(7V - 1)(7V - 2)^'^ + 4a + 2c + (77V - 12)^ + 4) (74) 

and in general 

_ + ac + (37V - 4)(7V - + (37V - 4)a/3 + 2(7V - l)c/3 + 4a + 2c + (77V -9)^ + 4 _ \ 

^ a + c+(27V-2)/3 + 3 ' 

n;+i (a+(iV-i)/3 + l) 

7 ^ T (75) 

n;^i (« + c+(2A^-i-l)/3 + 2) 

Moving further, 

(a + (AT - l)/3 + 1) /[O] P[3] 

^[3] = 7 sn TT ^ w w V 

(o + c + (2Ar - 4)/3 + 2j (a + c + (2Ar - 3)/3 + 2 j (a + c + (2Ar - 2)^ + 2 j (a + c + (2Ar - 2)/3 + sj (o + c + (2Ar - 2)/3 + 4j 
where 

P[3] = (a + 3)(a + 2)(a + c - /3 + 2)(a + c - 2/3 + 2)+ 
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+{N - l)/3((a + C + 2)(6a2 + 5ac + 32o + 11c + 40 + 4/3^) + 2/3(a + 2c + 2)(2a + c + 5)) + 

+(Af - l)(iV - 2)/32 (^(16a^ + 21ac + Sc^ + 80a + 52c + 98) + 2/3(20a + 16c + 51) + 36^^^ + 
+2{N -1){N -2){N -3)P^{10a + 2ip + 7c+26) + 10{N - 1){N - 2){N - 3){N - 4)l3^ 

• One can get rid of the non-factorizable polynomials in I[Y] by switching to peculiar linear combinations. 
For example, 

(N-1)B 1 + NB (a+(iV-l)/3+l)(a+(7V-l)/3 + 2) 

1 + P ^ + P (a + c+{2N -2)p + 2j(a + c+{2N -2)13 + 3) 



Note that not only a decomposition into Uncar factors is obtained in this way, also the factor + c + {2N — 

3)|8 + 2^, which was present in the denominators of both I [2] and /[II], is canceled in this combination. 
As noted in [45] the relevant linear combinations are actually the Jack polynomials P'-^^^^[Y]: 

n 

l<ii<i2<...<i„ fc=l 

Pi^m[2]=m,iz) + ^m^,.^{z)=^zf + ^ Y: z.zj, 

P i l<i<j 

P(V«[3] = ^j3j(,) + ^m[,,ii(.) + ^^_|^_^m[i3](.), p(i//^)[2, 1] = mp,i](.) + ^m(,3](z). 



...(78) 



Indeed, according to [60] , the Selberg integrals of Jack polynomials are factorized: 



and the (3- and AT-dependent coefficient is 

where h parameterizes the transposed Young diagram Y = {hi > ^2 >...}. In particular, 

^ NpjN^ + 1) ^ A^ 

(ni-l + ni/3)(ni-2 + n2/?) /3 + 1 ' ^ ^ 

in accordance with (77). 

The Selberg integrals I[Y] satisfy a set of sum rules. 

N 
i<3 



Since {I - Zi) = I - J2^=i + J2i'<j ^iZj - ■■■ one has 



J,+i[0] = /e[0] - Arj,[l] + ^ ^ > /c[ll] - ■ • • = E ^c[l"] (82) 

n=0 ' 



what is indeed true for (72). This sum rule involves only the single row Young diagrams. 

N 
i<j 

Je+2[0] = /c[0] - 27V/e[l] + iV/e[2] + N{N - l)7e[ll] - • • • (83) 



Similarly, from the expansion 11^(1 ^ ^i)'^ = 1 — 2 J^iLi Zi + Y^f=i zf + 2 Y^^<i ZiZj — ... one gets 



which includes only the double row Young diagrams (of which the single row diagram is a particular case with 

fc2 = 0). 

Similarly, expanding nr(-'^ ~ Zi)"^, one can deduce the expansion of /c+to[0] into a sum of the m'-iow Young 
diagrams with m' < m. Moreover, such sum rules can also be written for /c+m[5^] with arbitrary Y. 
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